***************************************
 * Descriptive Statistics *
***************************************

use "replication_data.dta"

* Summarise variables *
estpost summarize oecd anz_bc anz_activityoutlook nationalgovt cabinet_coalition govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu

* Export to RTF file *
esttab using "all_descriptive.rtf", cells("count mean(fmt(%4.3f)) sd(fmt(%4.3f)) min max") replace title({Descriptive Statistics.}\label{tab1})

* Export to TEX file
esttab using "all_descriptive.tex", cells("count mean(fmt(%4.3f)) sd(fmt(%4.3f)) min max") replace title({Descriptive Statistics.}\label{tab1})

***************************************
 * Trends Over Time *
***************************************

* Line Graphs for all Variables for Figure *

line oecd edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(96(2)104) title(OECD Business Confidence Index, size(vlarge)) ytitle("") xtitle("")
graph export time_oecd.jpg, replace

line anz_bc edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(-100(50)100) title(ANZ Business Confidence, size(vlarge)) ytitle("") xtitle("")
graph export time_anzbc.jpg, replace

line cabinet_coalition edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(0(0.5)1) title(Cabinet Coalition, size(vlarge)) ytitle("") xtitle("") 
graph export time_cabinetcoalition.jpg, replace

line anz_activityoutlook edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(-50(25)50) title(ANZ Activity Outlook, size(vlarge)) ytitle("") xtitle("") 
graph export time_anzao.jpg, replace

line annualgdpgrowthe11yearout edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(-5(2.5)5) title(Annual GDP growth 1 year out, size(vlarge)) ytitle("") xtitle("") 
graph export time_gdpgrowth.jpg, replace

line unemploymentrateg11yearout edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry)) yline(0, lcolor(white) lstyle(solid)) ylabel(0(2.5)12.5) title(Unemployment rate 1 year out, size(vlarge)) ytitle("") xtitle("") 
graph export time_unemployment.jpg, replace

line usdollarexchangerateh21yearout edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry)) yline(1, lcolor(white) lstyle(solid))  ylabel(0.25(.25)1) title(US dollar exchange rate 1 year out, size(vlarge)) ytitle("") xtitle("") 
graph export time_usd.jpg, replace

line australiandollarexchangeratei21y edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  yline(1, lcolor(white) lstyle(solid))  ylabel(.6(.2)1) title(Australian dollar exchange rate 1 year out, size(vlarge)) ytitle("") xtitle("") 
graph export time_aud.jpg, replace

line exports_billion edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(0(5)25) title(Exports of goods and services (nominal $ billion), size(vlarge)) ytitle("") xtitle("") 
graph export time_exports.jpg, replace

line housesales_1000 edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(10(10)40) title(House sales (thousand), size(vlarge)) ytitle("") xtitle("") 
graph export time_houses.jpg, replace

line v58 edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(-80(40)80) title(Perception of Monetary Conditions (net) 1 year out, size(vlarge)) ytitle("") xtitle("") 
graph export time_monetary.jpg, replace

line gepu edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))   ylabel(0(100)400) title(Global Economic Policy Uncertainty, size(vlarge)) ytitle("") xtitle("") 
graph export time_gepu.jpg, replace






/*

ssc install multiline, replace

multiline oecd nationalgovt annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion   edate, missing   xline(`=mdy(12,1,1990)' `=mdy(12,1,1999)' `=mdy(12,1,2008)' `=mdy(12,1,2017)' `=mdy(12,1,2023)')

* OECD Business Confidence *
twoway line oecd edate, yaxis(1) xline(`=mdy(12,1,1990)' `=mdy(12,1,1999)' `=mdy(12,1,2008)' `=mdy(12,1,2017)' `=mdy(12,1,2023)')|| line annualgdpgrowthe11yearout edate,  yaxis(2)  legend(size(vsmall))

generate upper = 104
local barcall upper y if inrange(y, 1990, 1999) | inrange(y, 2008, 2017) | inrange(y, 2023, 2024), base(96)
twoway bar `barcall' || line oecd edate yaxis(1) xline(`=mdy(12,1,1990)' `=mdy(12,1,1999)' `=mdy(12,1,2008)' `=mdy(12,1,2017)' `=mdy(12,1,2023)' lcolor(black))|| line annualgdpgrowthe11yearout edate, yaxis(2)

line oecd edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(96(2)104) title(OECD Business Confidence Index) ytitle("") xtitle("")

line anz_bc edate, xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)', lcolor(blue)) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)', lcolor(cranberry))  ylabel(-100(50)100) title(ANZ Business Confidence) ytitle("") xtitle("")


 `=mdy(12,1,2008)' `=mdy(12,1,2017)' `=mdy(12,1,2023)')|| line annualgdpgrowthe11yearout edate, yaxis(2)

line anz_bc edate if year>1998, yaxis(1) xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)')|| line annualgdpgrowthe11yearout edate if year>1998, yaxis(2) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)')  legend(size(vsmall))

line anz_bc edate if year>1998, yaxis(1) xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)')|| line unemploymentrateg11yearout edate if year>1998, yaxis(2) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)')  legend(size(vsmall))

line anz_bc edate if year>1998, yaxis(1) xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)')|| line gepu edate if year>1998, yaxis(2) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)')  legend(size(vsmall))

 line anz_bc edate if year>1998, yaxis(1) xline(`=mdy(12,1,1990)' `=mdy(12,1,2008)'  `=mdy(12,1,2023)')|| line australiandollarexchangeratei21y edate if year>1998, yaxis(2) xline(`=mdy(12,1,1999)'  `=mdy(12,1,2017)')  legend(size(vsmall))

line anz_bc edate, yaxis(1)|| line annualgdpgrowthe11yearout edate, yaxis(2)

line anz_bc annualgdpgrowthe11yearout edate

line anz_activityoutlook edate

set obs 1000
gen x = rnormal(0,1)
gen y = rnormal(0,2)
graph twoway (sc x y) || (scatteri 2 1 -2 1, recast(line))*/

***************************************
 * OLS Regression *
***************************************

* OECD 1987-2024 *

reg oecd nationalgovt 

estimates store oecd1 


reg oecd nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year

* reg oecd nationalgovt govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion

estimates store oecd2

coefplot (oecd1, keep(nationalgovt) mcolor(blue) ciopts(color(blue))) (oecd1, keep (cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion mcolor(black) ciopts(color(black)))),  xline(0, lpattern(dash) lcolor(black)) 

* 


* OECD 1999-2024 *

reg oecd nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu year

estimates store oecd3

coefplot (oecd2, keep(nationalgovt) mcolor(blue) ciopts(color(blue))) (oecd2, keep (cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu mcolor(black) ciopts(color(black)))),  xline(0, lpattern(dash) lcolor(black)) 


* ANZ BC 2006-2024 *

reg anz_bc nationalgovt 

estimates store anz1

reg anz_bc nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu year

estimates store anz2

coefplot (anz2, keep(nationalgovt) mcolor(blue) ciopts(color(blue))) (anz2, keep (cabinet_coalition govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu mcolor(black) ciopts(color(black)))) ,  xline(0, lpattern(dash) lcolor(black)) 

* Marginal Effects of Government & Coalition *

reg anz_bc nationalgovt#cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu 
margins nationalgovt#cabinet_coalition
marginsplot, level(90) xlabel(0 "Labour-led Government" 1 "National-led Government") xtitle("") recast(scatter)

reg anz_bc nationalgovt#cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu year
margins nationalgovt#cabinet_coalition
marginsplot, level(90) xlabel(0 "Labour-led Government" 1 "National-led Government") xtitle("") recast(scatter) ciopts(lwidth(medthick))

/*
reg anz_activityoutlook nationalgovt annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu if year<2020
estimates store pre2020

reg anz_activityoutlook nationalgovt annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu if year>2019
estimates store post2020

coefplot pre2020 anzoutlook,  xline(0, lpattern(dash) lcolor(black)) 
*/

/*Building the Table (RTF file) */
esttab oecd1 oecd2 oecd3 anz1 anz2 using "tablea1.rtf", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001) r2 title({Determinants of Business Confidence.}\label{tab1})

/*Building the Table (TEX file) */
esttab oecd1 oecd2 oecd3 anz1 anz2 using "tablea1.tex", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001) r2 title({Determinants of Business Confidence.}\label{tab1})


***************************************
* VIOLIN PLOTS *
***************************************

ssc install violinplot, replace
ssc install dstat, replace
ssc install moremata, replace
ssc install palettes, replace
ssc install colrspace, replace
ssc install grstyle, replace

gen oecdpos = oecd/100

gen oecd90 = oecd - 90

gen anzpos = anz_bc + 100


sum anz_bc
sum anzpos

* https://github.com/benjann/violinplot

ssc install vioplot
vioplot oecd,  over(nationalgovt) 


violinplot oecd, pdf(ll(0)) noline fill(line) over(nationalgovt) overlay legend(label(1 "Labour-led Government") label(2 "National-led Government"))


violinplot oecd90, pdf(ll(0)) noline fill(line) over(nationalgovt) overlay legend(label(1 "Labour-led Government") label(2 "National-led Government"))



violinplot anzpos, pdf(ll(0)) noline fill(line) over(nationalgovt) overlay legend(label(1 "Labour-led Government") label(2 "National-led Government"))


violinplot oecd , pdf(ll(0)) split(nationalgovt) key(fill)  

violinplot oecd, over(nationalgovt) overlay left dscale(.)  rag(offset(-.005 -0.01) mlcolor(%40))  

violinplot anzpos, over(nationalgovt) overlay left dscale(.)  rag(offset(-.005 -0.01) mlcolor(%40))  
	
violinplot anzpos , pdf(ll(0)) split(nationalgovt) key(fill)    

***************************************
* TESTS FOR AUTOCORRELATION *
***************************************

eststo clear

gen qdate = qofd(edate)
tsset qdate

* OECD 1987-2024 *

newey oecd nationalgovt, lag(1)
estimates store oecd1 

newey oecd nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion, lag(1)
estimates store oecd2


* OECD 1999-2024 *

newey oecd nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu, lag(1)
estimates store oecd3 

* newey oecd nationalgovt govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu, lag(0) 
* estimates store oecd4 




* ANZ BC 2006-2024 *

newey anz_bc nationalgovt, lag(1)
estimates store anz1

newey anz_bc nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu, lag(1)
estimates store anz2

* newey anz_bc nationalgovt govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu, lag(0)
* estimates store anz2


/*Building the Table (RTF file) */
esttab oecd1 oecd2 oecd3 anz1 anz2 using "rtable1.rtf", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001)  title({Determinants of Business Confidence.}\label{tab1})

/*Building the Table (TEX file) */
esttab oecd1 oecd2 oecd3 anz1 anz2 using "rtable1.tex", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001)  title({Determinants of Business Confidence.}\label{tab1})



***************************************
* TABLE A1 in Appendix *
***************************************

* OECD 1987-2024 *

reg oecd nationalgovt govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year
estimates store oecd1 


* OECD 1999-2024 *

reg oecd nationalgovt govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu
estimates store oecd2




* ANZ BC 2006-2024 *

reg anz_bc nationalgovt govt_seats annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu
estimates store anz


/*Building the Table (RTF file) */
esttab oecd1 oecd2 anz using "tablea1.rtf", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001) r2 title({Determinants of Business Confidence with Alternative Measure of Government Constraints.}\label{tab1})

/*Building the Table (TEX file) */
esttab oecd1 oecd2 anz  using "tablea1.tex", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001) r2 title({Determinants of Business Confidence with Alternative Measure of Government Constraints.}\label{tab1})
***************************************
 * Dominance Analysis *
***************************************


capture ssc install tuples

* OECD 1987-2024 *
domin oecd nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year, reg(reg) fitstat(e(r2))

* OECD 1999-2024 *
domin oecd nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58  gepu, reg(reg) fitstat(e(r2))

* ANZ BC 2006-2024 *
domin anz_bc nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58  gepu , reg(reg) fitstat(e(r2))

***************************************
 * ACTIVITY OUTLOOK ANALYSIS *
***************************************

* ANZ ACTIVITY OUTLOOK 2006-2024 *

reg anz_activityoutlook nationalgovt 
estimates store anzoutlook1

reg anz_activityoutlook nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu
estimates store anzoutlook2

/*Building the Table (TEX file) */
esttab anzoutlook1 anzoutlook2 using "outlook.tex", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001) r2 title({Determinants of ANZ Activity Outlook, 2006-2024.}\label{tab1})

coefplot (anzoutlook1, keep(nationalgovt) mcolor(blue) ciopts(color(blue))) (anzoutlook2, keep (nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu mcolor(black) ciopts(color(black)))) ,  xline(0, lpattern(dash) lcolor(black)) 

* ANZ OUTLOOK 2006-2024 *
domin anz_activityoutlook nationalgovt cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58  gepu , reg(reg) fitstat(e(r2))


***************************************
 * PM ANALYSIS INSTEAD OF PARTY *
***************************************

***************************************
 * OLS Regression *
***************************************

* OECD 1987-2024 *

reg oecd i.PM cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year
margins i.PM
marginsplot, xlabel(1 "Lange" 2 "Palmer" 3 "Moore" 4 "Bolger" 5 "Shipley" 6 "Clark" 7 "Key" 8 "English" 9 "Ardern" 10 "Hipkins" 11 "Luxon") recast(scatter) yline(100, lpattern(shortdash) lcolor(black)) title("OECD Business Confidence Index, 1987-2024", size(large)) xtitle("") ciopts(lwidth(medthick)) ytitle("Linear Prediction", size(medium))

estimates store oecd1 

coefplot, drop(cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu constant)  xline(0, lpattern(dash) lcolor(black)) 


* OECD 1999-2024 *

reg oecd i.PM cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu 
margins i.PM
marginsplot, xlabel(5 "Shipley" 6 "Clark" 7 "Key" 8 "English" 9 "Ardern" 10 "Hipkins" 11 "Luxon") recast(scatter) yline(100, lpattern(shortdash) lcolor(black)) title("OECD Business Confidence Index, 1999-2024", size(large)) xtitle("") ciopts(lwidth(medthick)) ytitle("Linear Prediction", size(medium))
estimates store oecd2 

coefplot, drop(cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu constant)  xline(0, lpattern(dash) lcolor(black))  


* ANZ BC 2006-2024 *

reg anz_bc i.PM cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion year housesales_1000 v58 gepu
margins i.PM
marginsplot, xlabel(6 "Clark" 7 "Key" 8 "English" 9 "Ardern" 10 "Hipkins" 11 "Luxon") recast(scatter) yline(0, lpattern(shortdash) lcolor(black)) title("ANZ Business Confidence, 2006-2024", size(large)) xtitle("") ciopts(lwidth(medthick)) ytitle("Linear Prediction", size(medium))
estimates store anz

coefplot, drop(cabinet_coalition annualgdpgrowthe11yearout unemploymentrateg11yearout usdollarexchangerateh21yearout australiandollarexchangeratei21y exports_billion housesales_1000 v58 gepu constant)  xline(0, lpattern(dash) lcolor(black)) 


/*Building the Table (RTF file) */
esttab oecd1 oecd2 anz using "pm.rtf", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001) r2 title({Determinants of Business Confidence.}\label{tab1})

/*Building the Table (TEX file) */
esttab oecd1 oecd2 anz using "pm.tex", b(3) replace unstack nobaselevels label nogap se star (+ 0.1 * 0.05 ** 0.01 *** 0.001) r2 title({Determinants of Business Confidence.}\label{tab1})
